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Outline 
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Nanotube Comparison 


Carbon Nanotubes *3?^ 

jjfc, ■*' j 

*<f**+l- j 

Boron Nitride Nanotubes -‘iv? 

Mechanical Properties 
(Axial Young’s Modulus) 

1 .33 TPa 
(very stiff) 

1.18 TPa 
(very stiff) 

Axial Thermal Conductivity 

>3000 W/mK 
(highly conductive) 

-300-3000 W/mK 
(highly conductive) 

Axial Coefficient of 
Thermal Expansion 

-1 xIO 6 K 1 
(very low) 

-1 x 1 0' 6 K' 1 
(very low) 

Electric Properties 

Metallic or semiconducting 

Insulator (6.0 eV band gap, 10 15 ' 18 
ohm-cm) 

Dielectric strength 2x1 0 5 V/mm 
Piezoelectric (0.25-0.4 C/m 2 ) 

Thermal Oxidation Resistance 

Stable up to 
300-400 °C in air 

Stable to over 
800 °C in air 

Neutron Absorption 
Cross-Section 

C = 0.0035 barn 

B = 767 barn (B 10 -3800 barn) 

N = 1 .9 barn 

(Excellent neutron radiation shielding) 

Color 

Black 

White 

(can be dyed to color) 
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Langley All-Nanotubes Actuator/Sensor Film 



CNT film electrode 


BNNT/Polymer active layer 
CNT film electrode 



Goal: flexible, mechanically durable, large actuation, high sensitivity 
Modeling is essential in achieving these goals 
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Molecular Dynamics Model of BNNTs 


• Define forces between atoms using a 
given interatomic potential (energy) 

V - Eh ('•«)- W ('•«)} F, - -W'M 

Uj 

• Evolve atoms according to Newton’s 

law: a i = F i /m i 



Tersoff- Brenner type of potential: 

1. Albeetal., (Rad. Eff. Def. Sol. 141 (1997)85) 

2. Verma - Sekkal: (Nanotechnology 18 (2007) 435711) 

Total potential energy: 
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where: 

D n 


Vr ( r ij) = F c ( r u)jz J exp (-/^V25 (r u -/■„)) - repulsive term 

n rv 

Va ( r v) = F c ( r ij)j^ ex p(-p'lys ( r ij - r o)) - attractive term 
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F ‘ W "| 0 - 5 1 + cos 
0 

B u = { l+ r np ;) 


( n r ij -R l \ 

TV— 

R 2 ~R\ 


r < R t 

, R l < r ; < R 2 - cut-off function 


^2 < r „ 


1 

2 n 


P ij = ^ F M k )G(6 ijk )exv\tf(r ij -r i ^ 


k*i,j 




d » d o+(K~ c ose m ) 


- angular dependent 3 -body term 


0 jjk is the angle between bonds i - j and i - k. 


MD models do not currently address the piezoelectric behavior 
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Tension 


Shear 


Piezoelectric Properties of BNNTs 

Piezoelectric Effect 


h-BN plane BNNT 



(Sai & Mele, PRB 68, (2003) 241405) 

Induced polarization, /?, under strain e jk : p i = e ijk e jk 
where e ijk - piezoelectric tensor with symmetry: 
e =-e =-e =-e ; e =0.086-0.12 e I Bohr 

xxx xyy yxy yyx ’ in 


The MD model has to reproduce the characteristic polarization behavior 



MD Model for a Piezoelectric BNNT 


Introduce a dipole term to the interatomic potential 


v-v,*u r -, u r - 2 

i={B}j={N} 


Dipole Parametric Equation: 




o 



p Q J d - fitting parameters 




N 


h-BN unit cell 


The dipole term is simple, efficient, and allows for analytic expression of the 

piezoelectric tensor 
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Analytical Model for a Piezoelectric BNNT 

Shell Theory of the Hexagonal BN Lattice* 


During homogeneous deformation the h-BN 
lattice sustains force equilibrium by shifting the 
N-sublattice at a shift vector £ with respect to the 
B-sublattice. 

5 is determined by the equilibrium condition: 

00/0 = 0 with 

® = 'z M v (r i] ;cos 6 m ,k*i,j) and 

%;cos 

as a function of the interatomic potential parameters: 


Shift between B and N sub-lattices 



8r r 


2 3 

with A = 1 — 


/ d 2 y ^ 

V^/o 


+ 12r n 


a 2 y 

ydr^d cos Ojjf, /o 


4n; 




and for r tj and cos 6 ijk : 


dr: 


+ 12r n 


y /o 


d 2 V 


ydryd cos 6 ijk /q 


\ / 

+ 18 


d 2 V 


dcosd ijk dcosd ijk /o 


\ / 

-9 


d 2 V 


dcosd ijk dcosd ifl /o 


\ / 

+ 12 


av 


5cos ^W 0 


^ = r o 2 (<V + 2£ „/i)( n « + Cd(«/i+^)- r o 4 (««+?«) (^ + ^ )] 2 Z 1 2 


cos »* -rr("2 +?■)(<*+?,)>< ^ + 2 ^ +^*„a»[ 3 (”I ♦eJK + £,)- 2 ('*ii +£„)(< +£,)-2(< +£)] 


VI 


*Song et al., Journal of Nanoscience and Nanotechnology, 8, 3774-3780, 2008. 
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Expressing Polarization* 


Parametric equation: p B N,ij=Po 


r BN,ij r 0 


+ ^ 2(2 +COS ^ fe ) 


Substituting r tj , and cos 6 ijk , and limiting to linear terms of £ ap , and r 0 K ap , 
excluding rotation, e = e yx : 



General equation expressing polarization 
of any deformed h-BN layer through the 
interatomic potential (A and r 0 ), and two 
fitting parameters, p 0 and f 0 . 


! . _L"_L J 



I 1-1 h 





mm 


*Yamakov et al., Comp. Mat. Sci., 95, 362-370, 2014. 
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Polarization of a hexagonal BN Plane 


In plane deformation: K ap = 0 

P- c = ^=^-[(l + ^)-3r 0 (l- A)]^ + A)-3r 0 (l- A)]e ; 

P y = “ 2 ^/f r 2 [(1 + - (1 “ ^)] £x y 

Piezoelectric tensor: 

Px = e xxx £ XX + ^/xv 

n = p p = p p 

* >’ yxy xy yyx yx 



Stretched h-BN Plane 


e xxx e xyy e y xy e y yx ^ [(^ + ^ ^)] 

Preserves the hexagonal symmetry 


Model Parameters | 



U r ] 

' P 



r o 

(nm) 

A 

E 

(TPa) 

P o 

(e/nm) 

h 

£xxx 

(e/Bohr) 

Albe et al. 

0.1462 

0.276 

1.11 

0.029 

-1.4 

0.0894 

Verma-Sekkal 

0.1457 

0.241 

1.48 

0.041 

-0.75 

0.0870 



£ 


xx 


N. Sai, EJ. Mele, PRB 68 (2003) 241405 
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Polarization of a BNNT 

Zig-zag (n 0) tube: K xl = 0; K . = 1 / R; r 0 2 / R 2 « 1 


Folded along x-axis 



R = nr 0 \[?>/2 ji 



MD simulation of a 
BNNT under stretch 


P = Px = iSrZ K 1 + ' 3t ° (1 " A) ] (l + + 64SR 2 (1 " 27 ‘ e - ’ 

Valid in the limit r 0 2 K 2 yy = r 0 2 / R 2 « 1 

No polarization in twist, s xy Spontaneous polarization ~ 1/R 2 

p sp = const / R 2 


MD Results 

Piezoelectric Polarization in Stretch 





Polarization of a BNNT 

Armchair tube: (n n) k xx = 1 / R; K yy = 0; r 0 2 / R 2 « 1 


Folded along y-axis 




P Py 2 + 


No polarization in stretch, £ yy 
No spontaneous polarization 


MD Results 

Piezoelectric Polarization in Twist 




Folded at angle cp 



Polarization of a BN Nanotube 

m\P> 


Chiral tube ( n m)with chiral angle, tan cp = 


m + 2n 


g n=^^[(l^)-3^(l-A)](l + v) 

P x = e n £ s +e u e t ZVJr ° 

e,4= '^7^[ (1+ ' 4) " 3 ' e(1 "' 4) ] 


MD Results 



(n m) 


Polarization in Stretch 


Polarization in Twist 



Summary 


Parametric equation: 


PBN,ij Po 


r BN , ij r 0 



Analytic expression of the piezoelectric coefficients through the interatomic potential 

= - V = ~e„ = - V = ^j[( ! + A ) “ (* _ A )] 

The model reproduces the piezoelectric properties predicted by first-principle 
The model predicts 1 / R 2 spontaneous polarization of BNNTs 


Limitations of the analytical model: 

Valid for small curvature (bending): tubes with r 0 2 /R 2 « 1, or R > 1 nm (most experimental 
BNNTs) 
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